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Introduction 

Convected  large  scale  structures  have  been  identified  in  many  different 
flows  such  as  wall  boundary  layers,  wakes,  and  jets.  In  each  case ,  the 
presence  of  the  coherent  structures  can  have  a  significant  effect  on  the 
overall  flow  properties  such  as  pressure  distributions  and  mixing  rates. 

In  order  to  successfully  predict  the  overall  properties ,  the  dynamics  of 
the  large  scale  structures  must  be  understood  and  reasonably  modeled.  This 
requires  first  thaz  the  structures  must  be  identified  in  experimental  data 
and  it  is  the  purpose  of  this  paper  to  propose  an  improved  method  to  do  so. 

The  difficulties  involved  with  searching  for  a  convected  vortex  structure 
in  a  general  flew  fceld  have  been  treated  by  Lugt  .  In  a  previous  paper,  the 
present  authors  discussed  a  method  of  c\  --nt  structure  identification  which 
was  based  on  a  Fourier  transformation .  2..  technique  proved  to  be  very 

'useful  in  isolating  tne  coherent  vortex  structures  and  hence  illuminating 
the  dynamics  of  those  structures . 
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The  present  identification  technique  appear*  to  allow  the  identification 


of  smaller  and  less  well  defined  structures.  It  is  based  on  the  concept  of 
decomposing  the  flow  (specifically  the  velocity  field)  into  various  components 
of  rotation ,  translation ,  shear,  ere.  and  then  determine  the  magnitude  and 
importance  of  each  of  these  elements . 

The  Model 

The  flow  model  is  organized  in  terms  of  the  data  cell  illustrated  in 
Figure  1.  At  each  corner  of  the  data  cell  several  pieces  of  information  may 
be  available .  In  this  case  only  the  velocities  at  these  positions  are 
considered  since  the  search  is  for  a  velocity  structure .  However,  the  same 
technique  would  be  applicable  to  the  search,  for  structure  in  the  passive 
contaminants;  for  example  temperature  or  density. 

The  two  components  of  velocity  az  each  corner  of  the  data  cell,  u^, 
v.,  u ^  /  v2,  Uj,  *,*,/  ,  and  v  ,  are  knewn  by  experiment.  Decomposing  these 

velocities  into  their  component  structures  requires  a  model  of  those 
structures .  A  set  of  mutually  orthogonal  descriptors  is  shown  in  Figure  2a 
and  consists  cf  rotation ,  translation  in  each  of  two  orthogonal  directions , 
expansion  in  the  same  directions  and  shear.  The  mutually  orthogonal 
descriptors  are  chosen  because  they  are  physically  realistic .  Other 
orthogonal  descriptors  could  be  chosen,  as  shown  in  Figure  2b.  They  are  not 
included  in  the  present  technique  because  the  present  six  descriptors  appear 
to  be  sufficient  from,  a  physical  point  of  view.  It  will  be  seen  in  the 
results  that  the  six  descriptors  are  adequate  to  identify  the  vortex  structure . 
Attempts  are  currently  underway  to  incorporate  all  eight  descriptors  and 
identify  any  potential  improvement  in  the  technique . 

The  model  aquation ,  in  vector  form,  is 

v  =  x  s’  +  e" 


(1) 


where 


V  =  velocity  at  the  corners  of  the  data  cell 

X  =  matrix  of  computational  weights 

8  =  parameters  of  rotation ,  translation,  etc. 

t.  =  error  between  the  model  and  the  measured  velocities . 


In  expanded  form  the  equation  is, 


The  velocities  are  measured  experimentally  and  are  the  data  to  be  analyzed 
for  coherent  structure .  The  6  vector  consists  of  the  modeled  parameters  shown 
in  Figure  2.  The  matrix  is  composed  of  the  computational  weights  which  relate 
the  flow  model  (Figure  2)  to  the  actual  measurements .  The  difference  between 
the  results  of  the  model  and  the  actual  velocity  measurements  is  attributed 
to  an  error  represented  by  the  e  vector.  Of  course,  if  the  model  had  eight 
coefficients  instead  of  six,  the  result  would  be  eight  equations  and  eight 
unknowns.  This  would  allow  no  room  for  error.  However,  this  would  only  model 
the  data,  not  the  process  because  the  actual  data  does  have  errors  and  only 
allows  six  physically  realistic  parameters  as  shown.  The  two  additional 
degrees  of  freedom  are  incorporated  into  the  non- zero  error. 

The  matrix  of  computational  weights  is  developed  by  observing  the  effect 
of  each  of  the  descriptors  on  the  velocities  in  the  data  cell.  The  lengths  of 


of  the  sides  of  the  rectangular  data  cell  are  A  and  3  and,  for  example,  the 
velocity  u ,  produced  by  the  rotational  component,  3^  is  proportional  to  the 
length  of  the  side,  A  and  is  in  the  -X  direction,  hence  the  negative  entry  -A 
in  the  first  matrix  entry.  Likewise  the  effect  due  tn  rotation  on  the  v 
component  is  proportional  to  the  side  B  and  is  also  negative,  hence  the  entry 
-B  in  the  first  column,  second  row.  The  effect  on  u due  to  the  translation 
in  the  same  direction  3^  directly  proportional  and  therefore  the  entry  1 
in  column  two,  first  row. 

It  should  be  noted  that  only  the  proportionality  is  of  importance .  The 
exact  relationship  between  the  descriptors  and  the  velocities  is  not  critical 
because  that  scale  can  be  (and  is)  absorbed  in  the  3  values  yet  to  be  found 
and  that  the  objective  here  is  to  identify  the  spacial  distributions  of  the 
various  descriptors  in  the  hope  that  the  distributions  will  clarify  the  flow 
structure.  That  is,  the  objective  here  is  to  develop  contours  in  space  which 
describe  the  various  2  terms.  It  is  these  cc.. tours  which  will  define  the 
coherent  structures  and  the  numerical  values  associated  with  those  contours 
are  not  as  important  as  the  shape  of  the  contours . 

The  Solution 

The  presents  tier,  of  the  solution  will  be  in  the  form  of  spacial  contour 
maps  which  will  then  allow  the  identification  and  location  of  coherent 
structures .  The  method  of  solution  is  based  upon  the  fact  that  the  optimum 
values  of  the  3  parameters  (i.e.,  those  identified  as  the  best  fit  of  the  model 
to  the  data)  occur  when  the  error  vector  e  is  of  minimum  length.  This,  in 
turn,  requires  that  the  error  vector  be  orthogonal  to  the  surface  containing 
all  of  the  vectors  X '  which  compose  the  matrix  of  computational  weights  in 
Eqn. (2) .  Therefore, ^ 

Xj-e  =  0  and  x'^'S  = 


0 


(3) 


or  in  general 


X  e  -  0  (4) 

But  the  error  is  the  difference  the  actual  measurements  V  and  those  which 
would  make  the  error  zero  V,  so 

X  (V-V)  =  0  (5) 

A  - - 

but  V  =  X&  ( since  e  =  0)  (6) 


so 

or 


X  V  -  X  XS  =  0 

A  ^ 

6  =  (x'x)  x'v 


where  the  6x6  matrix,  because  the  components  are  independent,  yields 


v'y  = 
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(7) 

(8) 


(9) 


.2  .2 


T f  A  +B  =  1  (which  presents  no  precise;  since  only  relative  values  of 
3  are  of  importance)  the  inverse  yields 


(X'X) 


-1 


=  k  I 


(10) 


so 


8=4-  x'v 


(ll) 


The  matrix  X  and  the  velocity  measurements  v  are  known  and  B  can  be 
easily  evaluated.  Results  for  both  synthetic  flows  and  actual  flows  are 
presented  in  the  following  sections. 

Synthetic  Flow  Results 

In  order  to  examine  the  contour  plots  of  the  orthogonal  descriptors  due 
to  idealized  flow  fields,  a  synthetic  vortex  was  created  and  subjected  to  the 
manipulation  of  Egn.  (11).  The  velocity  distribution  in  the  synthetic  vortex 
is  shown  in  Figure  3,  where  the  velocity  vectors  are  indicated . 

By  the  methods  described  above,  the  contours  for  fixed  values  of  B  are 

n 

shown  in  Figures  4-8  for  various  orthogonal  descriptors .  Figure  4  presents 
the  contour  map  for  the  rotational  component,  $  .  As  might  be  expected ,  most 
of  the  contour  is  rotationally  symmetric  and  the  most  rapid  changes  in  B  values 
occur  where  the  velocity  changes  are  most  rapid  as  seen  in  Figure  3.  The 
deviations  from  rotational  symmetry  occur  in  regions  of  little  change  of 
magmizude  and  are  caused  by  the  face  char  although  she  field  is  rotationally 
symmetric,  as  zhe  rectangular  daza  cell  is  moved  through  the  field,  the  values 
it  computes  are  symmetric  about  the  center  and  the  rectangular  axes  but  cannot 
he  onf-roji)  r~-rstior.ally  symmetric .  T'-is  is  due  to  finite  difference  representa¬ 
tion  of  a  rotational  field  by  a  rectangular  grid. 

The  only  place  in  the  contour  plot  that  presents  a  problem  for  the 
rectangular  grid  in  that  region  ;;here  the  velocitu  chances  are  small  so 
contour  lines  are  farther  apart.  In  a  real  flow  such  a  region  would  be  unlikely 
and  in  fact  would  also  be  of  little  interest  from  the  point  of  view  of  coherent 
structure .  Thus  the  square  grid  presents  no  problem  in  terms  of  flow  structure 
identification  since  coherent  structures  are  invariably  identified  with  large 
changes  in  flow  properties ,  at  least  relative  to  their  own  scale. 


The  translation  components,  8^  and  8^,  produce  rather  interes ting 

contours  as  shown  i ,,  Figures  5  &  6  respectively .  In  each  case  a  zero  value 

contour  passes  through  the  synthetic  vortex  center  and  parallel  to  the 

direction  of  the  translational  velocity  component  in  question .  This  is 

really  a  reflection  of  the  fact  that  the  vortex  is  perfectly  symmetrical 

since  it  is  synthetic  and  that  the  velocity  field  is  fixed  and  not  translating . 

The  contour  for  the  expansion  component  in  the  horizontal  direction ,  8^, 

are  shown  in  Figure  7  and  result  in  a  four  lobe  configuration.  The  contour 

plot  for  expansion  in  the  vertical  direction  is  rotated  by  90°  from  that  in 

Figure  7  and  thus  results  in  exactly  the  same  plot.  As  in  the  previous 

Figures,  the  presence  of  the  vortex  leads  to  a  very  distinctive  pattern. 

The  shear  component ,  8  ,  is  plotted  in  Figure  8  and  results  m  a  contour 

6 

map  very  simi lar  to  that  observed  for  the  extension  components  out  rotated  by 
45°  about  the  center  of  the  synthetic  vortex .  This  is  cor.sistant  with  the 
fact  that  the  maximum  would  occur  in  a  plane  inclined  at  45c  to  the  principal 
stress  directions. 

The  couuczz  plots  described  above,  due  to  the  synthetic  vortex,  will  be 
the  guide  for  the  structure  which  can  be  identified  in  real  flows  based  on 
the  ability  tc  recognize  the  appropriate  patterns . 

Real  Flow  Results 

a.  Vortex  Structure  near  a  Wall 

An  instantaneous  flow  field  containing  multiple  vortices  is  described  in 

Reference  4  and  is  produced  by  a  cam  shaped  rotor  which  produces  a  vortex  each 

time  it  rotates  at  the  boundary  of  the  flow  field.  The  study  of  the  vortex 

structure  is  motivated  by  some  attractive  potential  applications  of  the  rotor 
5-7 


geometry . 


An  example  of  the  structure  of  the  flow  field  produced  by  this  method  is 
shown  in  Figure  9.  This  field  has  been  analyzed  by  the  technique  of  Reference 
2  to  determi  ri'->  the  vortex  locations .  In  a  frame  of  reference  moving  with  the 
centers  <^r  the  location  of  the  two  vortices  found ,  the  coherent  structures 
be  seen  in  Figures  10  and  11,  showing  the  upstream  and  downstream  vortex 
structures  respectively .  It  is  the  objective  of  the  present  method  to  identify 
the  vortices  more  accurately ,  thereby  allowing  the  use  of  the  new  technique  to 
isolate  even  smaller  structures  in  complicated  flow  fields. 

The  rotational  orthogonal  descriptor ,  3^  of  Figure  2,  isolates  the  field 
of  Figure  9  and  is  shown  in  Figure  12.  According  to  the  results  of  the  synthetic 
vortex  rotational  component  shown  in  Figure  4,  the  presence  of  a  vortex 
structure  is  indicated  by  a  closed  contour  with  a  maximum  or  minimum  at  its 
center  (depending  on  the  sign  of  rotation) .  In  this  case  the  rotational 
component  clearly  identifies  two  vortices .  upstream  and  downstream.  Each  of 
the  structures  has  a  negative  sign,  indicating  a  cicckrwise  rotation .  The 
location  of  eaon  vortex  so  ider.tifi  ad-  coincides  with  the  locations  found  by 
flow  vi snail tazcur.  and  by  a  method  of  trial  and  error  which  consists  of  sub¬ 
tracting  various  translational  veiocutias.  unzul  an  apparent  vortex  structure 
is  evident. 

In  each  case  the  locations  of  vortices  are  more  accurate  than  those  found 
by  the  technique  of  Reference  2.  In  addition  the  improved  potential  of  the 
present  technique  is  also  evident  in  Figure  12  because  not  only  the  two  major 
vortex  structures  cm  isolated  but  an  additional  vortex  structure  rotating  in 
the  opposite  direction  is  evident  in  the  lower  center  of  the  Figure.  This 
vortex  structure  could  not  be  identified  by  the  previous  technique  described 
in  Reference  2. 

The  fact  that  this  small  cortex  structure  actually  exists  is  verified  by 
the  smoke  flow  visualization  photograph  of  Figure  13,  which  shows  only 


approximately  half  of  the  field  in  Fiqure  12.  The  smoke  streak  lines  verify 

the  presence  of  the  structure  and  also  indicate  the  strength  of  the  technique . 

The  origin  of  the  structure  itself  is  believed  to  be  due  to  the  production  at 

the  rotor  or  to  breeding  by  the  larger  vortex  structure  itself.  It  has  been 

8,9 

shewn  by  Francis,  et.  al .  that  under  certain  conditions  a  plate  oscillating 

perpendicularly  out  of  a  wall  can  produce  two  vortices  of  the  same  sign  for 

each  cycle  of  the  oscillation .  ^  similar  effect  was  found  for  the  flow 

10 

produced  by  the  rotor  of  different  velocity  and  length  scales  .  In  this  case, 
however,  the  smaller  vortex  is  of  opposite  sign  to  the  major  vortex  structure . 
Thus  the  breeding  of  a  secondary  vortex  by  the  primary  structure  is  a 
probability . 

Looking  back  to  the  orthogonal  descriptors  of  the  synthetic  vortex,  the 
pattern  produced  by  the  stream  wise  translational  descriptor ,  Figure  5,  shows 
a  line  in  the  translational  direction  through  the  ■ vortex  center  with  closed 
contours  above  and  below  that  position .  Such  a  pattern  is  observed  twice  in 
the  results  of  Figure  14;  one  is  at  the  location  of  the  upstream  vortex, 
although  only  the  upper  lobe  of  the  closed  contours  is  evident.  The  other  is 
at  the  location  of  the  downstream  vortex  where,  conversely ,  only  the  lower 
closed  lobe  is  found. 

It  is  of  interest  to  note  that  the  actual  contour  values  in  the  case  of 
the  synthetic  structure  yield  zero  values  on  the  line  passing  through  the 
vortex  center.  When  a  vortex  is  translating  through  a  flow  field,  a  trans¬ 
lational  flow  is  superimposed  on  the  vortex  structure  so  that  the  rortex 
center  lies  essentially  on  an  average  value  contour  if  the  flow  field  without 
the  vortex  is  approximately  uniform.  If  the  flow  field,  without  the  vortex, 
is  still  rapidly  varying  then  this  need  not  be  the  case. 


Translation  in  the  orthogonal  direction  can  also  be  employed  to  locate 
vortex  centers  as  shown  in  Figure  25.  As  in  the  case  of  the  synthetic  vertex 
(Figure  6)  the  structure  is  a  vertical  line  passing  through  the  vortex  center 
with  the  value  of  the  contour  roughly  equal  to  the  average  value  if  the  field 
without  the  vortex  is  relatively  uniform.  As  in  Figure  14,  the  downstream  vortex 
is  not  found  as  easily  due  to  the  fact  it  is  more  diffuse  and  less  concentrated. 
There  is  a  possibility  that  additional  information  may  be  extracted  from  the 
field  in  terms  of  the  direction  of  translation  of  the  structure .  Further 
analysis  is  proceeding  along  these  lines. 

For  both  the  stream  wise  and  transverse  expansion  descriptors ,  Figures 
16  and  17,  the  pattern  for  the  vortex  location  should  be  similar  to  that  produced 
in  Figure  7  for  the  synthetic  vortex.  Three-fourths  of  this  structure  is  found 
in  the  upstream  vortex  of  the  contour  plots  of  Figures  16  and  17.  The  ndssing 
lobs  is  the  upstream  one  farthest  from  the  wail,  although  there  it  some  evidence 
of  it  in  Figure  16  and  more  indication  in  ~imure  17.  The  reason  for  this 
apparently  missing  piece  of  contour  is  probably  due  to  the  flow  over  the 
vortex-  structure,  due  to  the  fact  that  the  structure  is  moving  at  a  velocity 
less  than  that  of  the  free  stream  velocitry  far  from  the  wall. 

The  structure  is  extremely  weak  for  the  downstream  vortex  and  suggests 
the  possibility  for  using  this  technique  to  quantify  the  strength  of  the  vortex, 
'lore  test  cases  will  be  required  to  calibrate  the  technique  and  t,,e  effect  of 
a  non  uniform  flow  on  which  the  vortex  flow  is  convected. 

In  a  similar  manner,  the  shear  flow  descriptor  also  adequately  describes 
the  upstream  vortex  location  but  not  the  downstream  vortex.  The  position  of 
the  lobes  is  rotated  by  45°  so  they  no  longer  lie  in  the  quadrants  but  rather 
on  the  axes.  Again,  three  out  of  four  lobes  are  present  in  the  upstream  vortex 
but  the  downstream  vortex  is  not  evident. 


Thus  in  the  real  flow  case  the  most  useful  descriptor  appears  to  be  the 
rotational  component.  Figure  12.  The  rotational  contour  is  capable  of  locating 
not  only  the  concentrated  upstream  vortex  but  also  the  weaker  downstream  vortex. 
The  usefulness  of  the  other  descriptors  (in  addition  to  isolating  the  rotational 
descriptor)  is  yet  to  be  determined . 
b.  Vortex  Structure  in  a  Jet 

The  typical  instantaneous  velocity  field  of  an  oscillating  jet  is  shown 
in  Figure  19.  These  flow  fields  have  been  shown  to  contain  convected  large 
scale  structures ,  an  example  of  which  is  shown  in  Figure  20.  The  present 
objective  is  to  isolate  those  structures  more  precisely  and  in  a  manner  which 
will  allow  the  identification  of  smaller  structures . 

Applying  the  same  orthogonal  descriptor  technique  discussed  in  earlier 
sections  results  in  the  identification  of  two  vortex  structures  as  shown  by 
the  contours  of  the  rotational  orthogonal  cescripzor  in  Figure  21.  The  upper 
and  larger  of  the  mwo  vortices  is  apprcxirareig  three  times  as  strong  as  the 
lower  vertex.  The  upper  vortex  has  been  found  by  other  methods,  i.e.  the 
method  of  Ref  erer.es  1  and  the  visual  teenni  cue  of  s-jbtractir.g  stream  wise 
velocity  and  observing  the  flow  structure  in  essentially  different  frames 
of  reference .  In  addition ,  the  vortex  structures  have  beer,  observed  in  the 
experimental  case  by  smoke  flow  visualizazior..  The  lower  weaker  vortex  has 
been  identified  by  none  of  these  techniques . 

The  stream  wise  translation  orthogonal  descriptor  is  shown  in  contour 
form  in  Figure  22.  As  was  the  case  in  the  synthetic  vortex  and  in  the  rotor 
flow  field  as  well,  the  translational  contours  appear  to  follow  the  general 
direction  of  convection  of  the  large  structure  itself.  The  transverse  trans¬ 
lational  component ,  Figure  23,  clearly  indicates  a  zero  transverse  velocity 
contour  which  is  easily  verified  by  the  velocity  vector  field  of  Figure  19. 


The  expansion  components  in  the  two  perpendicular  directions  are  shown 
in  Figures  24  and  25  and  do  not  yield  much  information  as  to  coherent  structure 
The  lack  of  identifiable  structure  is  probably  related  to  the  weakness  of  the 
vortices  relative  to  those  found  in  the  case  of  the  rotor  flow  field  previously 
described.  The  same  is  true  of  the  shear  orthogonal  descriptor  whose  contours 
are  shown  in  Figure  26. 

Expanded  Use  of  the  Present  Technique 

The  data  cases  discussed  in  the  present  treatment  all  involved  cpacial 
distributions  of  data.  One  could  also  extract  information  from  a  set  of 
fixed  probes  in  a  flow  field  (in  particular  a  boundary  layer  as  in  Figure  27)  . 
The  continuous  signals  generated  by  these  probes  with  time  are  shown 
schematically  in  Figure  28  for  two  com.por.enzs  of  velocity  from  each  of  the 
three  probes.  The  actual  velocity  change  with  time  will  depend  upon  the 
strength  and  locacicn  of  any  vortices  converted  throug  the  field.  With  a 
sufficient  number  cf  probes  (where  sufficient  is  not  yet  defined)  the  continu- 
ous  field  can  be  degenerated  into  a  grid  which  is  equivalent  to  those  in 
Figures  3,  9,  and  19. 

Thus  the  temporal  data  can  be  analyzed  in  the  same  fashion  and  by  the 
same  technique  as  the  date  treated  in  the  present  case.  This  would  then 
allow  the  determination  of  the  vortex  or  large  scale  structure  in  a  flow  field 
and  would  allow  the  identification  of  that  structure  by  a  small  number  of 
variables . 

Conclusions 

The  data  analysis  method  developed  in  Reference  2  was  successful  in 
the  identification  of  large  scale  coherent  flow  structures .  However, 
additional  accuracy  was  desired  from  the  point  of  view  of  a  tecnnique  which 
is  capable  of  isolating  small  structures  which  would  othenvise  escape 


detection.  Identification  of  these  smaller  structures  then  offers  the 
possibility  the  temporal  data,  as  opposed  to  spacial  data,  can  also  be  analyzed 
to  determine  the  vortex  structures  being  convected  in  a  boundary  layer.  The 
smaller  scale  structures  can  now  be  found.  The  temporal  data  case  is  currently 


being  examined. 
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Figure  2a.  Orthogonal  Descriptors  included 
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Figure  2b.  Orthogonal  Vesrriptors  not  included 
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Figure  3.  Velocity  distribution  in  the  synthetic  vortex. 
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Figure  10.  Structure  of  Upstream  Vortex  by  Technique  of  Ref. 
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Figure  11.  Structure  of  Downstream  Vortex  by  technique  of  Ref. 


Figure  12.  Rotational  Descriptor  for  the  Rotor  Flow  field 


GO  =  -f-  3000  rpm 


Figure  13.  Smoke  flow  visualization 
of  the  rotor  flow  field 


14.  streamwise  translation  descriptor  for  the  rotor  flow  field 


Stream  wise  expansion  descriptor  for  the  rotor  flow  field 


tor  for  the  rotor  flow  field 


Figure  19. 


Typical  instantaneous  flow  field  in  an  oscillating  jet. 


Figure  20.  Sn:oke  flow  visualization  of  a  vortex 
in  the  oscillating  jet  (over  c) 


Figure  23.  Transverse  component  of  the  translational  velocity 
orthogonal  descriptor  for  an  unsteady  jet  flowfi el, 


Figure  27.  Multiple  hot  wire  probes  in 
a  boundary  layer  containing 
a  vortex. 
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Figure  28.  Schematic  of  the  signals 
obtained  from  the  probes 
of  Figure  12. 
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Summary 

The  acoustic  field  produced  by  an  unsteady  jet  is  examined  both  in  the  free 
jet  configuration  and  whan  the  jet  is  placed  in  a  rectangular  duct.  Standing 
transverse  acoustic  waves  are  identified  in  the  duct  and  they  correspond  to 
transverse  duct  dimensions  which  lead  to  a  resonance  i^ith  the  driving  frequency. 
The  presence  of  the  strong  acoustic  waves  offers  the  possibility  of  employing 
these  waves  to  transfer  energy  by  means  of  normal  forces  as  opposed  to  the  more 
traditional  view  of  an  essentially  viscous  transfer  by  means  of  the  physical 
mixing  of  two  flows.  From  the  point  of  view  of  potential  application,  such 
driven  duct  flows  may  correspond  to  ejectors  or  to  jet  pumps.  Therefore,  the 
mass  flow  through  the  rectangular  duct  is  measured  and  analyzed  as  a  gross 
indication  of  the  momentum  transfer  from  the  jet  to  the  induced  flow. 

Introduction 

The  use  of  unsteady  flow  to  enhance  the  transfer  of  momentum  from 
a  primary  to  a  secondary  flow  has  been  pursued  by  various  investigators  operating 
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from  different  points  of  view.  Much  of  the  effort  has  centered  on  creating  a 

time  dependent  primary  flow  to  transfer  momentum  by  normal  (pressure)  forces 

rather  than  tangential  (shear)  forces.  Such  a  transfer  would  be  more  efficient 

and  could  be  accomplished  in  a  shorter  streamwise  distance.  A  secondary 

objective,  still  of  interest  even  if  the  primary  objective  is  not  attained, 

is  to  increase  the  mixing  rate  of  the  primary  jet  so  that  relatively  complete 

mixing  is  accomplished  in  a  short  streamwise  distance.  The  advantage  of 

1  2 

complete  mixing  for  ejector  applications  has  been  described  by  Quinn,  Bevilaqua 
and  others.  The  mixing  increase  must  be  achieved,  however,  without  an  exorbi¬ 
tant  increase  in  primary  nozzle  losses.  A  comprehensive  review  of  ejectors  has 

3 

recently  been  completed  by  Porter  and  Squyers  and  includes  unsteady  effects. 

The  objective  of  the  present  paper  is  to  outline  briefly  some  of  the  past 
efforts  in  unsteady  jets  and  duct  flows  and  to  highlight  some  current  work  in 
that  area. 

Some  Past  Unsteady  Experiments 

The  primary  objective  of  past  experiments  related  to  the  unsteady  transfer 

of  momentum  from  one  flow  to  another  has  been  to  increase  the  performance  of 

4 

ejector  devices.  An  early  experiment  was  performed  by  Lockwood  employing  a 
pulsejet  combustor  similar  to  those  used  in  German  aircraft  in  the  mid  1940's. 

The  pulse  consisted  of  a  slug  of  air  which  passed  through  the  shroud  in  the 
manner  of  a  piston.  An  analysis  of  such  a  system  where  the  primary  flow  (when 
it  is  on)  completely  fills  the  ejector  shroud  has  been  performed  by  Johnson 
and  Yang”*  by  the  method  of  characteristics.  These  results  are  in  reasonable 
agreement  with  the  results  of  an  accompanying  experiment  and  suggest  the  existence 
of  an  optimum  frequency. 

Although  the  concept  of  filling  the  duct  completely  with  a  slug  of  fluid 


is  enticing,  it  is  not  entirely  suitable  for  higher  area  ratios  between  the 


secondary  and  primary  flows.  For  such  cases,  Foa  proposed  a  rotary  primary 

jet  where  the  components  of  the  rotating  nozzle  were  pictured  to  act  like  vanes 

of  an  impeller.  The  rotary  jet  has  been  analyzed  by  Hohenemser  and  Porter^" 

8 

and  further  tested  for  application  by  Maise. 

Various  investigators  have  examined  different  ways  of  generating  a  time 

9 

dependency  in  a  free  jet  flow.  McCormack,  Cochran,  and  Crane  vibrated  the  jet 
nozzle  in  a  transverse  direction.  Binder  and  Favre-Marinet^  employed  a  rotating 
butterfly  valve  upstream  to  generate  a  pulsing  flow  at  the  exit  of  an  axisymmetric 
nozzle.  Flow  visualization  studies  by  Curtet  and  Girard^  verified  this 
structure  and  showed  that  the  pulsating  flow  resulted  in  a  series  of  vortex 

12 

rings.  This  pulsating  jet  was  subsequently  applied  in  an  ejector  configuration 

and  indicated  an  improvement  in  performance  especially  for  short  ejectors. 

Unfortunately,  the  performance  in  this  case  was  defined  as  the  Ejector  Thrust/ 

Nozzle  Thrust  with  a  fixed  total  pressure.  This  definition  does  not  take  into 

account  the  nozzle  thrust  efficiency  nor  any  changes  in  primary  mass  flow 

between  the  two  conditions,  so  the  actual  performance  is  still  in  question. 

A  fluidically  controlled  rectangular  nozzle  which  produces  a  flapping  jet 
13 

was  developed  by  V lets  and  exhibited  a  rapid  half  width  growth.  The  advan¬ 
tages  of  the  fluidic  control  are  the  lack  of  moving  parts  and  fluid  control 

14  15 

lines.  Detailed  investigations  of  this  jet  ~r*  showed  strong  vortex  inter¬ 
actions  leading  to  the  mixing  improvements.  Vortex  growth  in  driven  unsteady 
jets  at  low  Reynolds  Numbers  has  been  shown  by  Rockwell. ^  Another  nozzle 
employing  a  fluidic  instability  was  developed  by  Hill  and  Green^  and  consists 
of  a  rapid  expansion  section  in  a  round  duct. 


Simmons,  Platzer  and  Smith  have  examined  the  unsteady  jet  flowfield 


produced  by  an  angular  oscillation  of  the  jet  nozzle  itself  and  hav/e  recently 

19 

studied  a  jet  driven  by  an  embedded  oscillating  vane. 


In  addition  to  the  emphasis  on  the  primary  jet  itself,  some  recent  effort 

has  been  directed  at  the  control  of  wall  boundary  layers  by  unsteady  stimula- 
20  21 

tion  ’  which  has  application  to  the  mixing  d"ct  problem.  Again  the  aim  is 
to  create  vortex  structures  to  improve  mixing  and  energize  the  lower  part 
of  the  boundary  layer. 


Acoustj-C  Effects 

Host  of  the  unsteady  primary  nozzles  described  above  have  considered  the 

performance  advantage  of  such  nozzles  to  be  basically  due  to  the  improved  mixing 

rate  thereby  produced.  In  application,  this  requires  the  thrust  efficiencies 

of  the  unsteady  nozzles  to  be  high  enough  to  make  them  competitive  with  nozzles 
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such  as  the  hypermixing  type.  *  If  the  thrust  efficiencies  are  not  sufficiently 

high,  the  loss  in  the  primary  nozzle  cannot  be  overcome  by  a  viscous  mixing 

improvement.  There  is,  however,  an  alternate  phenomenon  which  can  be  employed 

to  improve  the  mixing  by  transferring  the  energy  through  normal  (pressure) 

forces  rather  than  tangential  (viscous)  forces  as  discussed  in  Section  1.  The 

mixing  can  be  accomplished  by  an  acoustic  interaction  between  the  primary  jet 
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and  the  duct  va_ls  as  demonstrated  by  Quinn.  His  configuration  consisted 

of  an  axisymmetric  duct  and  a  choked  primary  jet,  which  produces  a  strong 
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screech  tone  as  shown  by  Powell  and  recently  studied  by  Krothapalli, 
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Baganoff,  Hsia  and  Karamcheti. 

The  acoustic  interaction  is  most  effective  when  the  acoustic  frequency 

available  from  the  primary  jet  is  well  matched  with  the  fundamental  natural 
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modes  in  the  duct  shroud.  As  demonstrated  by  Viets,  Campbell  and  Korkan, 
a  substantial  difference  in  frequency  between  the  jet  and  the  fundamental 
modes  greatly  reduces  the  acoustic  effect  and  the  resultant  mixing  improvement. 
For  a  rectangular  duct  with  a  rectangular  primary  nozzle,  matching  the  acoustic 
frequencies  can  become  a  problem  if  one  simultaneously  desires  a  reasonable 


high  secondary  to  primary  flow  area  ratio.  In  addition,  for  jet  flows  of  less 
than  unit  mach  number  at  the  exit,  there  is  no  distinct  frequency  generated 
but  rather  a  broad  band  noise.  Therefore,  under  conditions  of  an  unchoken 
primary  jet  or  a  choked  primary  in  conjunction  with  a  larger  secondary  to 
primary  area  ratio,  an  alternate  acoustic  generator  is  required  if  an  acoustic 
interaction  is  to  be  provoked. 

Une  possibility  tor  use  as  an  unsteady  primary  jet  under  these  conditions 

is  the  "whistler  nozzle"  described  in  Reference  17,  since  it  can  operate  at 

13 

subsonic  mach  numbers.  The  acoustic  feedback  fluidic  jet  is  also  a  possi¬ 
bility  if  sufficiently  high  frequercy  can  be  generated.  A  rotor  jet  based  on 
the  mechanical  devices  described  in  Refs.  20  and  21  Is  -1jo  a  possible  choice 
for  a  primary  and  offers  the  additonal  advantage  of  allowing  a  number  of  primary 
nozzles  to  have  not  only  the  same  frequency  but  also  a  fixed  phase  relation¬ 
ship  and  is  currently  being  sf"died.  This  is  also  true  of  a  modified  version 

13 

of  the  acoustic  feedback  jet  nozzle  described  in  Ref.  26  as  a  gust  tunnel 
application. 

A  Subsonic  Acoustic  Generator 

Still  another  nozzle  which  offers  the  required  acoustic  properties  at 
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subsonic  exit  velocities  was  originally  developed  as  a  two  phase  injector, 
operating  with  a  liquid  primary  into  a  gaseous  secondary.  The  nozzle  has 
the  advantage  of  no  moving  parts  and  hence  high  reliability.  Operating  in  a 
single  phase  condition  in  air,  the  nozzle  produces  distinct  frequencies  which 
can  be  tuned  both  by  changes  in  the  nozzle  geometry  and  by  the  operating 
conditions.  The  nozzle  schematic  is  shown  in  Figure  1.  The  nozzle  operation 
simply  requires  that  the  jet  attach  to  one  wall  within  the  nozzle  body.  Then 
part  of  the  jet  is  scooped  off,  fed  back  and  impressed  on  the  jet  flow  just 
downstream  of  the  contraction.  This  causes  the  jet  to  attach  to  the  opposite 


wall  and  the  process  begins  anew.  The  strongly  unsteady  quasi-sinusoidal 
flowfield  is  clear  in  Figure  2  which  illustrates  the  two  phase  operation  of 
water  flow  into  air. 

The  remainder  of  this  paper  will  consider  the  flowfield  generated  by 
this  nozzle  and  its  acoustic  interaction  with  the  duct  geometry.  If  the 
acoustic  interaction  is  substantial,  then  the  nozzle  can  be  optimized  in 
terms  of  scale,  frequency  and  spreading  race  relative  to  a  particular  duct 

geometry  to  allow  its  application  to  an  ejector-type  device. 
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The  acoustic  interactions  observed  by  Quinn  *  were  attributed  to 
standing  transverse  waves  in  the  duct  due  to  a  tuning  of  the  screech  tones 
produced  by  the  choken  jet  to  the  transverse  dimension.  This  interaction 
probably  causes  the  jet  itself  to  oscillate.  In  the  present  case,  the  jet 
is  driven  unsteady  even  without  an  acoustic  interaction.  Thus  the  effect  of 
the  acoustics  on  the  jet  may  be  rather  small. 

Free  Jet  Acs us  tics 

Three  unsteady  jets  and  one  slot  jet  were  investigated  in  this  program. 

Two  of  the  unsteady  nozzle  designs  are  shown  in  Figure  3.  The  primary 
difference  between  the  two  designs  is  the  length  of  the  feedback  loop.  The 
different  designs  lead  to  a  frequency  difference  of  roughly  a  factor  of  three 
at  nozzle  stagnation  pressure  ratios  of  about  1.15  based  on  ambient  pressure. 
Nozzle  3  is  essentially  the  same  design  as  Nozzle  2  except  for  a  somewhat  shorter 
feedback  loop  which  increases  the  frequency'  by  about  20%. 

In  order  to  correlate  the  flow  inside  the  oscillating  jet  with  the 
acoustic  structure  outside,  hot  wire  measurements  within  the  jet  were  compared 
with  pressure  transducer  measurements  and  microphone  measurements  taken  outside 
the  jet.  A  Flow  Corp.  (now  Dataraetrics)  Model  900  constant  temperature  hot 


wire  anemometer  was  employed  inside  the  jet  and  an  Endevco  Model  8506-2  high 
frequency  response  pressure  transducer  as  well  as  a  B&K  Model  4136  microphone 
and  amplifier  system  were  employed  outside  the  jet.  The  output  from  each  of 
the  systems  was  input  to  a  Hewlett-Packard  Model  3582A  spectrum  analyzer  which 
yields  the  magnitude  as  a  function  of  frequency.  Thus  the  dominant  frequencies 
present  in  any  signal  could  be  easily  determined.  A  typical  spectral  output 
for  Nozzle  #2  at  a  stagnation  pressure  ratio  of  1.17  is  shown  in  Figure  4  as 
determined  with  the  Endevco  pressure  transducers. 

The  two  lowest  frequency  modes  present  in  the  flowfield  on  Nozzle  //2  are 
compared  in  Figure  5  as  they  are  determined  both  in  the  jet  (hot  wire)  and  out 
of  the  jet  (microphone  and  pressure  transducer).  The  apparent  frequency  steps 
in  the  hot  wire  case  are  only  caused  by  the  finite  bandwidth  employed  by  the 
spectrum  analyzer.  The  conditions  inside  and  outside  the  jet  flow  clearly 
indicate  the  same  frequency  content.  The  pressure  is  non-dinensionalized  with 
the  ambient  pressure  while  the  frequency  is  reduced  to  a  Strouhal  number  form 
by  V/D  (exit  velocity/exit  dimension) . 

The  dominant  frequencies  outside  the  jet  flow  are  compared  in  Figure  6 
for  the  microphone  and  pressure  transducer  results.  In  each  case  the  instrument 
is  located  in  -he  exit  plane  of  the  jet  and  at  a  distance  equal  to  half  the  major 
dimension  of  the  ejector  shroud  (in  other  words,  located  at  the  same  position 
as  the  shroud  wall,  when  the  shroud  is  mounted) . 

Each  of  the  frecuency  modes  A-G  can  be  identified  and  explained.  Mo',Q  A 
is  the  driving  mode  of  the  jet  oscillation.  Mode  B  is  the  second  harmonic, 
being  double  the  frequency  of  Mode  A.  Mode  C  is  very  weak.  In  some  of  the  spectra 
it  hardly  appears.  The  fourth  harmonic  is  Mode  E,  while  Modes  F  and  G  are  again 
only  multiples  of  Modes  E  and  F  respectively.  The  intermediate  Mode  D  can  be 
explained  based  on  the  nozzle  results  in  3311  Hz,  which  is  the  Mode  D.  Since 


the  Mode  D,  based  on  the  nozzle  scale  is  so  close  to  Mode  C,  the  third 
harmonic  of  the  driving  motion,  the  energy  from  Mode  C  appears  to  be  con¬ 
tributing  to  Mode  D  and  thereby  reducing  its  own  amplitude. 

The  frequency  of  Nozzle  //3  is  increased  over  that  of  Nozzle  #2  as  shown 
by  Figure  7.  The  first  harmonic  exists  at  approximately  1300  Hz  or  a  Stroudal 
number  of  .059.  At  higher  pressures,  the  second  harmonic  clearly  exists  at 
approximately  2600  Hz  (i.e.  curve  B)  and  the  third  harmonic  at  curve  D.  The 
intermediate  mode,  curve  C  is  caused  by  the  duct  mode  based  on  the  depth  of 
the  rectangular  nozzle,  5.08  cm.  as  was  the  case  in  Nozzle  #2. 

As  the  pressure  is  lowered,  the  frequency  of  Mode  B  rises  while  the 
frequency  of  the  duct  Mode  C  decreases.  The  two  modes  coalesce  by  P  =  1.08. 

At  higher  pressures  a  third  harmonic  (Mode  D)  is  evident  but  weak,  while  at 
lower  pressures  a  higher  frequency  value  appears.  The  fourth  harmonic  exists 
at  low  pressures  while  at  higher  values  Mode  E  appears  to  be  the  multiple  of 
the  average  of  Modes  B  &  C. 

Free  Jet  Flowfields 

As  night  be  ejected  based  on  previous  flapping  free  jet  experiments 
(Refs.  13,  15,  13  &  27),  the  unsteady  jet  mixes  more  rapidly  with  the  ambient 
fluid  than  jets  without  a  time  dependent  driver.  The  centerline  velocity  decay 
and  the  half  width  growth  of  Nozzles  1  &  2  are  shown  in  Figures  8  and  9 
respectively  and  are  compared  to  the  two  dimensional  slot  jet  value.  Based 
on  the  velocity  decay.  Figure  8,  there  appears  to  be  a  small  frequency  depend¬ 
ence  showing  improved  mixing  (decay)  as  the  frequency  decreases.  This  result 
is  opposed  to  that  of  Reference  19  which  found  a  mixing  improvement  with 
f  requency . 

The  two  apparently  conflicting  results  can  fit  into  a  coherent  picture 
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of  the  vortex  growth  in  the  unsteady  jet.  As  clearly  shown  by  Brown  '  and 


others,  a  forced  unsteady  jet  leads  to  large  vortex  formations.  For  the  case 
of  the  flapping  jet,  these  vortices  have  been  found  by  both  flow  visualization 
and  by  conditioned  sampling  with  a  hot  wire  anemometer.^  For  low  frequencies, 
these  vortices  would  be  widely  spaced  and  weak  and  exert  little  influence  on  the 
the  entrainment.  Under  these  conditions,  an  increase  in  frequency  would  improve 
the  mixing  properties.  On  the  other  hand,  at  higher  frequencies  there  is  little 
to  be  gained  by  having  the  vortices  even  more  closely  spaced. 

The  effect  of  vortex  growth  may  also  be  seen  in  the  halt  width  growth 
results  of  Figure  9.  For  each  frequency,  there  appears  to  be  a  plateau  or 
leveling  of  the  half  width  growth  at  a  different  streamwise  location.  The 
most  obvious  case  is  the  frequency  of  980  Hz  and  occurs  at  approximately  13 
nozzle  width  downs  t  =>m .  For  higher  frequencies,  this  occurs  progressively 
farther  downstream  and  for  f  =  3500  Hz  the  plateau  is  twice  as  far  downstream 
as  for  980  Hz.  Apparently  the  near  field  entrainment  depends  upon  the  frequency. 
Brown  has  observed  that  the  vortices  entrain  in  the  near  field  but  do  not  entrain 
external  fluid  farther  downstream.  Thus  the  plateaus  may  reflect  the  end  of 
the  near  field  entrainment  caused  by  the  vortex  structure. 


Duct  Acoustics 

A  rectangular  duct  was  chosen  as  the  ejector  shroud  as  shown  in  the 
schematic  of  Figure  10.  The  inside  dimensions  of  the  duct  are  5.1  cm  by  17.8 
cm.  The  duct  length  was  varied  during  the  experiment.  The  longer  transverse 
dimension  of  the  duct  was  chosen  to  achieve  resonance  with  a  driving  frequency 
of  approximately  1000  Hz. 
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As  described  by  Succi,  et  al.  and  employed  by  Viets  et  al.,  the 

frequencies  of  the  modes  corresponding  to  resonant  transverse  waves  are 

related  to  the  speed  of  sound  Cq  and  the  duct  dimensions  b  and  d  by 


F 

m,n 


C 

o 

2 


(m/b)2  +  (n/d)2 


where  m  and  n  are  integer  values. 


If  the  speed  of  sound  is  assumed  to  be  335  m/sec  (1100  ft/sec),  then  the 
frequencies  corresponding  to  some  particular  resonant  modes  are  shown  below  for 
duct  dimensions  5.1  cm  by  17.8  cm. 

m  n  F(Hz) 

1  0  946 

2  0  1887 

3  0  2833 

0  1  3289 

0  2  6578 

1  1  3422 

2  1  3791 


It  should  be  noted  that  these  values  do  not  take  into  account  the  flow 
in  the  duct  nor  any  real  fluid  or  real  wall  effects.  Actual  resonant 
frequencies  are  -nr; acted  to  approximate  these  values  but  are  not  precisely 
equal  to  them. 

In  order  to  test  the  duct  acoustics,  white  noise  (produced  by  a  General 
Radio  Type  1390  A  'Random  Noise  Generator)  was  introduced  at  the  duct  inlet 
and  the  resulting  accustic  levels  measured  at  the  exit  plane  of  the  duct 
(located  at  L/D  =  6.1).  A  typical  spectrum  of  the  acoustic  output  on  the 
duct  exit  centerline  is  shown  in  Figure  11a.  The  clearly  defined  peaks  indi¬ 
cate  a  response  of  the  duct  to  those  frequencies.  Calibration  of  the  level 
of  the  white  noise  source  shows  that  the  peaks  between  150  and  400  Hz  are 
exaggerated  out  of  nroportion  due  to  the  fact  that  the  input  levels  are  high. 
The  lowest  transverse  resonant  frequency  was  shown  to  be  approximately  946Hz. 
The  peaks  between  400  Hz  and  the  lowest  transverse  frequency  are  attributed 
to  the  organ  pipe  modes  described  by  Morse  and  Ingard, 
mode  of  150  Hz  (for  this  length). 


which  have  a  lowest 


If  the  acoustic  measurements  are  taken  off  the  duct  centerline,  even  to 
the  point  of  reaching  the  far  wall  location,  the  peaks  between  400  and  900  Hz 
are  not  greatly  affected  as  seen  in  Fig.  lib.  This  is  consistant  with  their 
dependence  on  the  duct  length.  However,  the  modes  near  the  1000  Hz  transverse 
resonance  value  do  vary  with  transverse  position  as  might  be  expected. 

To  test  the  longitudinal  modes,  the  end  of  the  duct  is  filled  with  sound 
absorbing  foam  which  should  eliminate  the  organ  pipe  mode  by  now  allowing 
reflections  of  the  waves.  That  this  is  accomplished  is  seen  in  Figure  11c 
(with  the  microphone  inserted  through  the  foam)  especially  in  comparison  with 
Figure  11a.  It  is  interesting  that  there  remains  a  small  peak  at  990  Hz 
(remember  that  the  apparent  low  frequency  modes  are  caused  by  the  input  signal. 

Transverse  Waves 

The  overall  objective  of  the  oscillatory  duct  flow  is  to  examine  the 
transverse  wave  structure  to  see  if  such  waves  are  actually  generated.  High 
frequency  response  Endevco  transducers  are  mounted  to  the  duct  walls  at  a 
streanwise  position  equal  to  two  duct  dimensions  (35.6  cm)  .  The  output  signals 
are  fed  to  a  Hewlett-Packard  Model  358A  Spectrum  Analyzer  which  decomposes 
the  signal  into  its  frequency  components. 

(a)  Resonance  in  a  long  duct,  X/D  =  6.1,  Nozzle  #2 

The  variation  of  each  component  magnitude  across  the  duct  is  shown  in 
Figure  12  for  Nozzle  # 2,  having  a  lowest  harmonic  frequency  approximately 
equal  to  the  anticipated  transverse  resonant  duct  frequency.  The  first 
harmonic,  1060  Hz,  produces  a  standing  wave  of  a  half  wavelength  with  a  node 
at  the  center.  The  second  harmonic  yields  an  entire  wave  with  two  node  points. 
The  third  and  fourth  harmonics  are  very  consistant  with  1%  and  two  waves 
respectively . 


The  same  procedure  applied  across  the  short  side  of  the  rectangular  duct 
yields  the  result  of  Figure  13  which  clearly  shows  that  the  standing  wave 
structure  is  not  a  function  of  that  direction. 

Taking  the  same  pressure  measurements  at  an  upstream  position,  x  =  .28 
(near  the  ejector  inlet),  and  decomposing  the  signals  into  their  frequency 
components  results  in  Figures  14  and  15.  The  structure  of  the  first,  second, 
and  fourth  modes  are  similar  to  the  downstream  structure  (Figures  12  and  13) 
except  for  several  data  points  near  the  center  of  the  duct.  This  discrepancy 
is  caused  by  the  primary  jet,  which  has  a  high  velocity  at  that  streamwise 
location  and  interferes  with  the  pressure  measurement.  An  additional  discrepancy 
between  the  upstream  and  downstream  modes  is  shown  in  the  two  third  (3200  Hz) 
modes.  The  problem  is  resolved  by  the  fact  that  3200  Hz  is  the  transverse  duct 
mode  corresponding  to  the  nozzle  depth  as  was  discussed  in  relation  to  the  free 
jet  inodes.  Figures  6  and  7. 

(b)  Resonance  in  long  duct,  X/D  =  6.1,  Nozzle  #1). 

The  spectra  developed  by  the  pressure  traverse  of  the  duct  at  S/D  = 

2  downstream  exhibit  a  strong  peak  at  the  primary  frequency  of  Nozzle  #1, 

3300  Hz,  as  shown  in  Figure  16.  Based  on  the  results  of  Figure  12,  a  waveform 
of  lh  waves  is  expected  at  the  3300  Hz  frequency.  However,  only  a  full  wave 
is  found.  There  may  be  an  effect  due  to  the  nozzle  duct  frequency  which  is 
almost  equal. 

The  next  substantial  wave  structure,  at  7800  Hz,  shows  no  definite 
pattern.  However,  if  the  resonant  frequency  of  the  longer  transverse  dimension 
of  the  duct  is  1060  Hz  (from  Figure  12) ,  a  corresponding  resonance  in  the  short 
duct  dimension  is  3700  Hz.  The  7800  Hz  frequency  is  roughly  double  that  value, 
so  one  might  expect  a  waveform  similar  to  that  of  the  2120  Hz  case  in  Figure  12, 
but  in  the  short  dimension.  That  is  indeed  a  possible  waveform  as  seen  in  Figure  17. 


At  frequencies  of  10,300  and  15,500  Hz,  higher  order  mode  shapes  are 
possible.  The  former  could,  for  example,  correspond  to  the  mode  M  =  5,  n  =  2 
while  the  latter  may  correspond  to  M  >  5,  n  =  2.  Of  course  as  the  order  of  the 
modes  increases,  it  is  more  difficult  to  determine  the  precise  modes  since  the 
spectrum  becomes  very  crowded  with  possible  modes.  From  the  point  of  view  of 
the  unsteady  ejector,  the  higher  order  modes  are  really  of  less  interest 
because  they  are  expected  to  contain  less  energy. 

(c)  Resonance  in  long  duct  X/D  =  6.1,  Nozzle  #3 

The  lowest  dominant  frequency  for  this  nozzle,  1300  Hz,  effectively 
produces  resonance  in  the  duct  according  to  Figure  18  and  verified  by  the 
uniform  magnitude  in  the  other  direction.  Figure  19. 

The  next  appreciable  wave  is  at  4640  Hz  and  is  not  a  uarLuiic  of  the 
first  mode  but  is  related  to  both  duct  dimensions.  From  Figures  18  and  19, 

M  =  2  and  n  =  1 .  The  calculated  frequency  for  this  mode  structure  is  3791  Hz, 
which  is,  of  course,  for  the  zero  flow  condition. 

The  6000  Hz  case  does  not  indicate  a  clear  wave  shape .  It  is  a  mixed 
node  depending  on  both  transverse  dinensions  and  is  almost  a  combination  of 
the  lower  frequencies.  It  may  be  a  result  of  an  oscillation  among  several 
modes  which  exist  near  this  frequency.  Additional  minor  mode  peaks  exist  at 
higher  frequencies  and  there  it  becomes  increasingly  difficult  to  distinguish 
between  modes. 

Thus,  when  the  forced  oscillation  is  near  the  resonant  frequency,  it 
achieves  a  manner  of  resonance  which  is  not  as  strong  as  it  would  be  if  the 
frequency  were  better  matched.  It  does,  however,  produce  similar  pressure 
waves  . 

(d)  Resonance  in  shorter  ducts 

The  standing  waves  described  above  appear  even  for  very  short  ejector 
shrouds.  The  transverse  waves  in  an  ejector  of  length  equal  to  its  major 
width,  i.e.  L/D  =  1,  arc  shown  in  Figure  20  driven  by  nozzle  1/2.  The  mode 


structure  is  essentially  the  same  as  in  Figure  12  and  shows  the  same  primary 
nozzle  interference  as  in  Figure  14 . 


Hass  Augmentation 

In  order  to  examine  the  potential  application  to  an  ejector  device,  the 
mass  augmentation  was  determined  by  integrating  the  exit  velocity  profiles  and 
dividing  by  the  primary  mass  flow  as  determined  by  the  venturi  measurement. 

The  mass  flow  ratios  are  shown  in  Figure  21  as  a  function  of  duct  length  for 
the  slot  nozzle  and  three  different  oscillating  nozzles.  The  actual  levels 
of  mass  augmentation  are  not  significant  because  the  nozzles  have  not  been 
optimized  in  terms  of  thrust  efficiency  and  spreading  rate  relative  to  the  duct 
size. 

The  trends  in  mass  flow  ratio  show  that  for  a  long  duct  the  use  of 
unsteady  flow  in  general  and  acoustic  interaction  in  particular  is  not  warranted. 
The  loss  in  nozzle  thrust  efficiency  suffered  by  the  unsteady  jets  has  over¬ 
powered  any  mixing  and  acoustic  advantages.  However,  for  shorter  duct  lengths, 
the  performance  does  not  deteriorate  as  rapidly  with  length  for  the  oscillating 
jets  as  compared  to  the  slot  jet.  It  is  for  these  very  short  ducts  that  the 
unsteady  effects  discussed  here  are  expected  to  be  of  significance. 

The  monotonic  decrease  in  performance  with  decreasing  length  is  only 
interrupted  at  a  length  of  2.1  shroud  widths.  This  position  will  be  examined 
in  greater  detail  but  the  preliminary  indications  are  that  there  is  no  strong 
streamwise  organ  type  of  resonance.  That  is,  there  is  a  coherent  streamwise 
wave  but  it  is  not  greater  in  magnitude  than  the  wave  which  exists  for  a  duct 
length  of  L/D  =  3. 


Conclusions 


The  fluidically  controlled,  positive  feedback  jet  has  exhibited  the 
potential  to  generate  a  significant  acoustic  interaction  between  parallel 


walls .  The  frequencies  can  be  controlled  by  the  nozzle  design  in  order  to 
.match  the  scale  of  the  duct  geometry.  The  resonant  condition  (i.e.  a  match 
between  the  jet  oscillation  frequency  and  the  natural  frequency  of  the  wall 
separation)  yields  a  strong  interaction  and  standing  transverse  waves  between 
the  walls.  The  existence  of  this  interaction  expected  to  be  a  prerequisite 
for  a  significant  energy  transfer  through  the  acoustic  mechanism.  Actual 
proof  of  this  transfer  will  require  additional  work. 

The  transverse  waves  generated  by  the  uns'eady  jet  have  potential 
application  to  ejectors,  both  for  thrust  augmentation  and  pumping.  The 
performance  of  such  devices  remains  to  be  optimized.  In  particular,  the 
design  of  the  nozzle  geometry  must  take  into  account  not  only  the  matching  of 
the  duct  geometry  to  the  frequency  but  also  simultaneously  matching  the 
spreading  rate  to  the  duct  geometry.  In  order  to  accomplish  this  result, 
preliminary  nozzle  experiments  are  required  to  determine  the  frequencies  and 
breeding  rates  of  a  number  of  various  designs.  From  these  a  design  can  be 
chosen  which  satisfied  both  the  acoustic  and  mixing  requirements. 
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Figure  2.  Two  phase  operation  of  Fozzle  1,  water  flow  into  ai: 
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figure  12.  Pressure  level  distributions  of  the  lower 
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Abstract 

The  three  dimensional  deformation  of  a  line  vortex  near  a  flat  wall  is 
examined  by  means  of  flow  visualization  and  quantitative  conditioned  sampling 
of  a  hot  wire  anemometer.  The  vortex  is  observed  to  undergo  a  very  rapid 
deformation  from  the  rectilinear  form  to  a  very  high  aspect  ratio  hairpin¬ 
like  shape.  Similar  shapes  have  been  observed  in  the  case  o^  driven  flows  in 
a  laminar  boundary  layer.  In  the  present  case  the  vortices  are  driven  by  a 
rotor  configuration. 

General  Background 

The  dynamics  of  vortex  structures  has  been  of  interest  to  many  researchers 
for  a  variety  of  reasons.  Much  of  the  early  interest  was  analytical  and 
motivated  by  the  ability  to  model  the  structure  in  terms  of  inviscid  hydro¬ 
dynamics. ^  The  ideal  ^ortex  structures  were  even  employed  at  one  point  as  a 

2 

model  for  atomic  structure  and  successrully  predicted  some  of  the  earlier 

observations  of  atomic  physics.  Multiple  vortex  interactions  wer  a  also 
3  4 

computed  ’  and  in  some  cases  stability  analyses  were  performed  to  reveal 
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potential  naturally  occuring  vortex  configurations,  for  example  the  series  of 
axisymmetric  ring  vortices  shown  in  Figure  1.  Of  course  the  most  successful 
application  of  inviscid  vortex  structures  to  predict  observed  a  .phenomenon  is 
the  well  known  Kavman  Vortex  Street  analysis."*  Karman  was  able  to  predict,  by 
actually  computing  the  effect  of  a  number  of  fortices  on  each  other,  that  only 
a  very  specific  spacing  of  vortices,  shown  in  Figure  2,  would  yield  a  stable 
configuration. 

Most  of  the  early  vortex  analyses  and  experiments  concerned  free  vortices 
(i.e.  those  removed  from  wall  boundaries).  This  allowed  the  structure  to  be 
treated  inviscidly  by  methods  outlined  by  Lamb'*'  and  others^.  Most  of  these 
vortex  structures  were  ring  vortices  and  these  have  been  reviewed  by 
Zaroodny ^ . 

The  dynamics  of  vortex  ring  motion  has  been  analyzed  inviscidly  by  Arms 

g 

and  Hama.  An  experimental  and  analytical  treatment  was  performed  by  Viets 

9  10 

and  Sforza  who  employed  a  viscous  model  due  to  Rankine  and  a  time  dependent 

11  12 

viscous  model  due  to  Hamel  and  Oseen  .  More  recent  fully  viscous  analyses 

,13  14  15 

are  due  to  Widnali  J  and  Saffman  and  have  been  reviewed  by  Saffman  and  Baker 
The  objective  of  the  present  paper  is  to  examine  the  three  dimensional 
vortex  dynamics  near  a  wall.  The  motivation  behind  looking  at  the  methods 
employed  for  the  prediction  of  vortex  dynamics  far  from  a  wall  is  simply  that 
these  methods  will  be  directly  applicable  '  .  is  willing  to  model  the  wall 

in  an  inviscid  fashion.  Thus  the  inviscid  wall  would  amount  to  a  plane  of 
symmetry  in  the  analysis. 

Vortex  Near  a  Wall 

The  initial  studies  of  vortex  structure  near  a  wall  were  driven  by  the 

interest  in  boundary  layer  transition.  The  initial  stability  predictions  of 
16  17 

Tollmien  and  later  Schlichtir.g  were  verified  by  the  classical  experiment 


of  Schubauer  and  Skramstad  who  employed  an  electromagnetically  oscillating 
ribbon  in  a  laminar  flow.  The  magnitude  and  frequency  of  the  oscillation  were 
controlled  as  well  as  the  Reynolds  number  of  the  flow.  The  results,  in  terms 
of  a  stability  diagram,  are  very  well  predicted  by  Tollmien's  analysis  (see 
Reference  16,  pg.  460). 

Even  more  interesting,  from  the  point  of  view  of  vortex  dynamics  in  the 

vicinity  of  a  wall,  is  the  structure  of  the  unsteady  flow  produced  near  the 

19 

surface.  As  calculated  by  Schlichting  for  a  neutral  (i.e.  neither  growing 

nor  decaying)  disturbance,  the  closed  streamlines  exhibit  a  vortex-like 

structure  which  is  convected  in  the  streamwise  direction  and  shown  in  Figure 

3.  This  structure  has  been  named  Tollmien-Schlichting  waves.  Some  excellent 

photographic  evidence  of  this  structure  behind  a  driven  oscillating  ribbon  of 

18 

the  type  employed  by  Schubauer  and  Skramstad  has  been  obtained  by  Wortmann, 
who  employed  ink  and  water  to  generate  the  streakline  pattern  shown  in  Figure  4 
The  deformation  process  of  a  straight  vortex  line  of  the  Tollmien- 
Schlichting  type  has  been  illustrated  by  ninze  ~  and  is  shown  in  Figure  5. 

The  vortex  distorts  slightly  and  this  small  distortion  is  amplified  by  the 
self-induced  velocity  distribution  set  up  by  the  vortex.  The  region  of  highest 

g 

curvature  of  the  vortex  has  the  highest  self-induced  velocity  .  Thus  the 
loop  in  the  vortex  is  propelled  fastest  from  the  wal]  and  thus  finds  itself 

in  a  higher  streamwise  velocity  region.  The  loop  is  then  convected  faster  in 

the  streamwise  direction  than  the  straight  portion  deeper  within  the  boundary 
layer  and  the  vortex  filiment  is  stretched  until,  by  a  mechanism  which  is  not 
fully  understood,  it  is  said  to  "burst".  This  mechanism  is  currently  being 
called  into  question  as  the  actual  process  leading  to  boundary  layer  transition 
Some  recent  evidence  indicates  that  the  initially  turbulent  region  (or  "spot") 

is  actually  composed  of  numerous  vortex  filiments  instead  of  a  single  deformed 


vortex.  However,  from  the  present  point  of  view  of  vortex  dynamics  near  a  wall, 
the  model  still  offers  a  rational  deformation  process  in  a  boundary  layer. 

Some  especially  good  photographs  of  the  deformation  of  an  initially 
rectilinear  vortex,  again  behind  a  vibrating  ribbon,  have  been  obtained  by 
Wortmann"*  employing  the  pulsed  hydrogen  bubble  technique  and  an  example  is 
shown  in  Figure  6,  a  view  of  the  deformed  vortex  loop  from  above  (looking  down 
perpendicularly  to  the  plate) . 

21 

In  the  deformation  model  proposed  by  Hinze  ,  shown  in  Figure  5,  the 

expected  streamwise  velocity  profiles  through  the  centerline  of  the  deformation 

are  shown  at  various  streamwise  positions  and  compared  to  the  typical  boundary 

layer  type  of  profile.  Actual  velocity  profiles  of  this  type  have  been  measured 
22 

by  Muiswinkel  ,  again  behind  an  oscillating  ribbon,  and  are  shown  in  Figure  7 
for  various  times  during  one  oscillation  of  the  ribbon. 

The  objective  of  the  portion  of  the  study  described  on  the  following 
pages  is  to  identify  the  deformation  of  an  initially  rectilinear  vortex  structure 
near  a  wall.  The  two  techniques  employed  in  the  study  are  flow  visualization 
and  conditioned  sampling  of  a  hot  wire  anemometer  signal. 

Flowfield  Geometry 

Various  techniques  are  possible  to  generate  the  initially  rectilinear 

18 

vortex.  The  vibrating  ribbon  employed  by  Schubauer  and  Skramstad  is  one 

possibility  which  has  been  discussed.  However,  in  the  present  experiment, 

the  objective  is  to  operate  at  a  larger  scale  than  the  previous  vibrating 

t  tbbon  experiments.  Another  poss.'bility  is  a  technique  employed  by  Francis 
23  24 

et.  al.  *  in  their  study  of  vorticity  generated  on  the  upper  surface  of 
an  airfoil.  The  basic  mechanism  is  a  plate,  positioned  perpendicularly  to 


the  surface,  which  alternately  extends  from  and  retracts  into  the  surface. 

Each  time  the  plate  is  extended  it  creates  a  recirculation  region  behind  itself 


and  this  region:’ is  shed  and  convected  downstream  when  the  plate  is  retracted 
into  the  surfaced 


From  the  present  point  of  view,  it  is  desirable  to  have  control  over  the 

strength  of  the  vorticity  produced  and  thus  the  flat  plate  is  rejected  in 

favor  of  the  rotor  geometry  shown  in  Figure  8.  The  rotor  has  discontinuity 

in  its  radius  at  the  cusp  shape.  The  mean  flow  is  from  the  left  and  the  rotor 

& 

turns  in  a  counter  clockwise  direction.  The  rotor  tip  appears  periodically 


above  the  surface  and  each  time  the  cusp  tip  is  exposed  to  the  mean  flow,  a 


ft 


transverse  vortex  (turning  in  the  clockwise  direction)  is  produced  lying 

•  "IS?' 


parallel  to  thej/j^urface  and  perpendicular  to  the  mean  flow.  The  strength  of 
the  vortex  so^raroduced  is  proportional  to  the  relative  velocity  between  the 
rotor  tip  and  the  mean  flow.  Thus  control  over  the  vortex  strength  can  be 


achieved  by  simply  changing  the  rotational  spaed  of  the  rotor,  including  the 
possibility  ofirtuming  the  rotor  in  a  clockwise  direction. 


The  rotor ^geometry  in  Figure  8  has  been  employed  to  study  the  vorticity 
"  25 

produced  as  a  tprasi— two-dimens ional  structure  .  Only  centerline  data  was 
taken  and  the  imaj ectory  of  the  vortices  produced  is  shown  in  Figure  9  for  a 


rotational  speed- of  w  =  3000  rpm  counter  clockwise.  Both  axes  have  been 
nondimens io naif zed  with  respect  to  the  maximum  height  that  the  rotor  extends 
above  the  surface,  H  in  Figure  8. 

In  addition  to  a  detailed  examination  of  the  vortex  structure  produced 

m 

by  the  rotor, ithe  geometry  has  been  applied  to  various  problems  as  shown  in 


.,26 


,  avoiding 


Figure  10,  including  the  control  of  the  flow  over  an  airfo 

27 

separation  over  a  rearward  facing  ramp  ,  driving  the  flow  in  a  dump  combustor 

-  &  '  29 

and  studying  the  properties  of  convected  structures  .  Each  of  the  potential 


applications  will  be  strongly  effected  by  the  three  dimensional  nature  of  the 

•  x 


vortex  flow 


Since  the  data  shown  in  Figure  9  was  obtained,  the  wind  tunnel  has  been 
rebuilt  and  usable  test  section  length  tripled.  The  objective  is  to  determine 
the  three  dimensional  character  of  the  vortex  structure  and  to  follow  that 
structure  farther  downstream  then  previously  possible. 

Flow  Visualization  Results 
a.  The  choice  of  rotor  speed. 

The  visualization  of  the  flowfield  is  accomplished  by  the  use  of  entrained 
kerosene  vapor.  The  kerosene  is  vaporized  by  slowly  dripping  it  on  an  inclined 
heated  flat  plate  where  it  is  vaporized  and  forced  through  a  series  of  tubes 
which  end  at  the  inlet  to  the  wind  tunnel.  Here  the  vapor  is  entrained  into 
the  tunnel,  yielding  the  streaklines  of  the  flowfield. 

The  initial  results  in  the  rebuilt  wind  bunnel  are  shown  in  Figures  11  - 
14,  at  the  same  conditions  as  tested  in  Reference  25,  a  rotational  speed  of 
w  =  3000  rpm  in  the  counterclockwise  direction. 

The  results  have  been  obtained  at  specific  phase  angles,  for  example 
Figure  11  shows  the  flowfield  for  a  phase  angle  0=0°  which  means  the  rotor 
is  in  the  upright  or  fully  extended  position,  as  may  be  seen  in  the  first 
photograph  of  the  Figure.  The  photographs  are  arranged  such  that  the  fields 
overlap  each  other.  Thus  portions  of  the  photographs  that  lie  under  each  other 
indicate  a  region  of  overlap. 

The  first  conclusion  drawn  from  Figures  11-14  is  that  the  scale  of  the 
vortex  disturbance  is  rather  large.  The  transverse  scale  of  the  tunnel,  from 
the  lower  wall  to  the  top,  is  30  cm.  In  Figure  11,  the  scale  of  the  vortex  at 
a  position  of  roughly  one  transverse  scale  downstream  can  be  estimated  to  be 
35  cm.  That  this  estimate  is  accurate  may  be  seen  from  Figure  10  of  Reference 
25.  Farther  downstream  in  Figure  10  (for  a  phase  angle  of  0°)  the  scale  becomes 
substantially  larger.  The  same  effect  may  be  seen  in  the  results  at  a  later 


phase  angle  (6  =  120°)  in  Figure  11.  The  newly  formed  vortex  near  the  rotor 
has  a  scale  of  perhaps  one-fourth  of  the  transverse  value  but  grows  very 
rapidly.  In  the  yet  later  phase  angles  the  development  of  the  vortex  may  be 
followed  in  the  streamwise  direction  and  it  is  clear  that  the  overlapping 
photographs  yield  the  same  structure  in  either  case. 

There  are  two  problems  with  the  3000  rpm  results  from  the  point  of  view 
of  studying  the  three  dimensional  vortex  structure.  One  is  that  the  scale  is 
too  large,  making  it  difficult  to  follow  the  details  of  the  rather  diffuse 
vortex.  The  second  is  that  the  vortices  are  relatively  closely  spaced,  being 
only  one  transverse  scale  apart.  This  might  create  some  problems  in  detecting 
influences  of  one  vortex  on  another.  For  these  reasons  it  was  determined  that 
a  lower  frequency  would  be  preferable.  In  addition,  a  well  formed  vortex  was 
desired  which  led  to  the  study  illustrated  in  Figure  15. 

At  a  phase  angle  of  0  =  30°  the  rotor  position  is  such  that  the  cusp 
point  is  rotated  30°  from  the  upright  position.  For  various  rotational  speeds, 
the  flowfields  are  shown  in  Figure  15,  varying  between  -600  and  -3000  rpm. 

In  each  case  it  is  clear  that  a  vortex  rotating  In  the  clockwise  direction  is 
produced,  indicating  that  the  freestream  velocity,  nominally  12  ra/sec,  is 
greater  than  the  tip  speed  of  the  cusp  shape.  A  possible  exception  is  the 
case  at  -3000  rpm  where  the  vortex  is  very  weak  or  non-existent  since  the  tit 
speed  is  approximately  equal  to  the  nominal  free  stream  velocity. 

From  Figure  15,  the  choice  of  a  rotational  speed  is  made.  The  objective 
is  a  v.’all  formed  vortex  with  a  scale  that  is  reasonable  for  the  tunnel  scale. 

A  rotation  speed  of  -1600  rpm  is  chosen  on  this  basis, 
b.  Vortex  generation. 

The  vortex  generation  process  employing  the  rotor  turning  in  the 
clockwise  direction  at  2000  rpm  and  the  flow  from  the  left  is  shown  in 
Figure  16  for  various  phase  angles .  At  zero  degrees  the  curling  up  of  the 


smoke  streaklines  behind  the  cusp  shape  is  clear.  A  secondary  frequency  is 
evident  on  the  smokeline  and  may  be  related  to  an  instability  of  the  flow. 

As  the  phase  angle  increases,  the  vortex  moves  downstream  and  at  a  phase  angle 
of  45°  the  vortex  interacts  with  the  corner  of  the  downstream  plate  behind 
the  rotor.  By  a  phase  angle  of  67.5°  the  vortex  is  more  clearly  defined  and 
by  90°  the  vortex  is  no  longer  evident.  This  is  the  first  indication  of  the 
three-dimensionality  of  the  flow,  since  the  plane  in  which  the  photographs  are 
taken  is  not  the  centerline  of  the  tunnel  but  rather  a  plane  closer  to  the  far 
wall. 

The  data  planes  considered  in  the  present  investigation  are  shown  in 
Figure  17.  In  addition  to  the  centerline  plane,  there  are  planes  at  five 
and  ten  centimeters  to  either  side  of  the  centerline.  The  photographs  of 
Figure  16  are  taken  in  plane  -4  (i.e.  four  inches  from  the  centerline).  The 
reason  that  the  vorcex  appears  to  disappear  into  a  cloud  of  smoke  at  a  phase 
angle  of  90°  is  simply  that  the  vortex  axis  has  shifted  to  an  orientation 
closer  to  perpendicular  to  the  wall. 

c.  Flowfield  fcr  1=0° 

Figures  13-22  illustrate  the  streaklines  in  the  flowfield  for  a  phase 
angle  of  0°  (i.e.  the  rotor  in  the  fully  extended  position)  and  in  the  various 
planes  shown  iu  Figure  17.  In  each  plane  the  curling  behind  the  rotor  is 
evident  and  in  the  centerline  case.  Figure  20,  a  vortex  is  clear  at  70  cm 
downstream  and  approximately  12  cm  off  the  wall  (See  vortex  A) .  Looking  for 
this  vortex  in  the  +  2  and  +  4  planes  shows  in  the  +  planes,  the  presence 
of  the  vortex  is  clear  while  in  the  +  4  planes  only  the  secondary  effect  of 
the  vortex  is  evident,  i.e.  the  curved  streaklines  in  Figure  18.  Thus  the 
vortex  is  perpendicular  tc  the  plate  in  the  +  2  planes. 

d.  Flowfield  for  0  =  90° 

For  a  phase  angle  of  90°  the  streaklines  in  various  planes  are  shown 
in  Figures  23-27.  Thr  vortex  located  immediately  behind  the  rotor  in  the 


last  sequence  has  now  moved  to  position  B  in  the  centerline  photograph  of 
Figure  25.  In  the  +  2  planes,  the  effect  of  this  vortex  may  be  seen  while 
in  the  +  4  planes  the  effect  is  small  but  present.  Thus  the  vortex  appears 
to  assume  a  three  dimensional  shape  very  rapidly,  in  this  case  within  a 
quarter  rotation  of  the  generator.  The  downstream  vortex  (”C")  has  moved  to 
80  cm.  downstream  and  15  cm.  above  the  wall.  Again,  the  effect  of  this  vortex 
is  very  small  in  the  +  4  planes  but  quite  clear  in  the  +  2  planes. 

e.  Flowfield  for  0  =  120°  and  0  =  180° 

The  vortex  presence  is  located  by  D  &  E  at  0  =  120°  and  F  at  0  =  180° 
in  Figures  28-37.  In  each  case  the  vortex  has  moved  farther  downstream,  so 
that  at  0  =  180°  only  one  vortex  is  evident  in  the  field,  approximately  30 
cm.  downstream.  The  presence  of  the  weak  vortex  lying  right  on  the  wall 
itself  is  clear  in  the  -4  plane  at  approximately  25  cm.  downstream.  Thus  the 
dynamics  of  the  vortex  process  involve  a  violent  curling  up  a^d  distortion  of 
the  control  portions  of  the  originally  rectilinear  vortex  but  a  weak  remainder 
of  the  vortex  which  continues  to  nave  near  the  wall. 

f.  Flowfield  for  8  =  240°  40=  270° 

By  a  phase  angle  of  240°,  the  vortex  in  the  centerline  plane  has  moved  to 

a  streamwi.se  position  of  40  cm  and  a  height  of  10  cm.  as  denoted  by  position 

in  Figures  38-42-  The  presence  of  the  vortex  is  strongly  felt  in  the  centerline 

plane  (Figure  40)  and  still  clear  in  the  +  2  planes.  Farther  from  the  center- 

line  only  a  slight  effect  is  seen,  consistent  with  a  weak  vortex  still  moving 

near  the  wall.  Some  evidence  of  a  secondary  vortex  ("H")  may  be  seen  at 

approximately  90  cm.  This  effect  may  be  due  to  the  fact  that  at  certain  speeds 

the  rotor  may  shed  more  than  one  vortex  per  cycle.  This  effett  was  observed 

during  some  studies  of  a  rotor  positioned  on  the  back  of  a  van  as  well  as  on 

the  back  of  a  car^.  The  secondary  vortex  has  also  been  seen  in  the  case  of 
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an  oscillating  orthogonal  plate 


The  same  general  effect  is  evident  in  the  case  of  a  phase  angle  of  270° 
Figures  43-47.  The  vortex  ("I")  is  identified  at  a  streamwise  position  of 
50  cm,  5  cm  above  the  plate.  Again  the  vortex  is  strong  on  the  centerline, 
smeared  at  the  +  2  positions  and  weak  and  near  the  wall  away  from  the 
centerline. 

g.  Conclusions  based  on  flow  visualization. 

The  flowfield  structure  depicted  by  the  visualization  results  of 
Figures  18-47  can  be  summarized  by  Figure  48.  The  initially  rectilinear 
vortex  deforms  very  rapidly  as  it  is  convected  downstream.  Initially  there 
is  some  evidence  that  the  vortex  lines  extend  out  to  the  front  and  back  walls 
of  the  tunnel.  However,  farther  downstream  the  connection  is  tenuous  and  it 
is  possible  that  the  vortlcity  lines  end  on  the  surface  of  the  flat  plate 
instead  of  the  wind  tunnel  walls. 

It  is  probably  that  the  single  deformed  vortex  is  due  to  the  end  conditions 
on  the  wind  tunnel  walls.  That  Is,  if  the  rotor  were  infinitely  long,  the 
initially  rectilinear  vortex  would  still  deform  into  a  similar  pattern  as  that 
in  Figure  48  and  that  the  pattern  would  be  repeated  with  a  certain  scale  length 
which  would  be  determined  by  a  question  of  the  stability  of  the  rectilinear 
vortex  to  an  applied  disturbance.  In  the  real  experimental  case,  the  applied 
disturbance  is  the  end  walls  which  supply  a  boundary  layer  which  will  cause 
the  vortex  ends  to  move  downstream  more  slowly  than  the  vortex  center.  This 
initial  deformation  then  progresses  into  the  shapes  observed.  The  aspect 
ratio  of  the  deformed  vortex  (i.e.  its  height  above  the  plane/distance  between 
the  legs  perpendicular  to  the  plate/  is  approximately  unity  by  30  generator 
heights  downstream. 
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Figure  3.  Streamline  pattern  associated  v/ith  a 
neutral  disturbance.  Calculations  by 
Schlichting1.9 
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Figure  7.  Instantaneous  velocity  profiles  at  various  times 
during  one  cycle,  at  one  streanwise  position, 
employing  an  oscillating  ribbon  (due  to 
Kuiswinkel22) . 
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Figure  9.  Trajectory  of  the  vortex  in  the  centerline 
plane,  w  =  3000  rpra  ccw 


Figure  13.  Centerline  flow  visualization  results  for  w  =  3000 
rpm  counterclockwise  and  phase  angle  0  =  150°. 
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Figure  15.  Formation  (at  0  =  30  )  of-the  vor£ices  formed  at  various  rotor  speed?, 
all  in  Che  clockwise  direction. 
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